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1 Introduction: 

The purpose of this note is to extend Dynkin isomorphim involving functionals 
of the occupation field of a symmetric Markov processes with non polar points 
and of the associated Gaussian field to a suitable class of non symmetric Markov 
processes. This was briefly proposed in [5] using Grassmann variables, extending 
to the non symmetric case some results of [3]. Here we propose an alternative 
approach, not relying on Grassmann variables that can be applied to the study 
of local times, in the spirit of [6]. It works in general on a finite space and 
on an infinite space under some assumption on the skew symmetric part of the 
generator which is checked on two examples. 

2 The finite case 

In this section, we will prove two formulas given in Proposition 1 and Corollary 2 
which relate the local time field of a non symmetric Markov process on a finite 
space to the square of the associated complex "twisted" Gaussian field. The 
result and the proof appear to be a direct extension of Dynkin's isomorphism. 

2.1 Dual processes 

Let us first consider the case of an irreducible Markov process on a finite space 
X, with finite lifetime generator L and potential V = (— L)~ x . 

Let m = [iV for any nonnegative probability /ionI. Recall that L can be 
written in the form L = q{I— II) with q positive and II submarkovian . Then the 
m-adjoint L can be expressed similarly with the same q and a possibly different 
submarkovian matrix II. Moreover, m = fiV, with jl the law of xq- under P M . 

Note that for any z = x+iy 6 C x , the "energy" < Lz,~z > m = X) ~{I j z) x ~z x m i 
is nonnegative as it can be written C x . y (z x — z y ){z^ — ~z^)+ < III + III — 
2, z~z > m ), with C x , y — C V:X — m x q x H x ^ y . The highest eigenvector of ^(11 + II) 
is nonnegative by the well known argument which shows that the module con- 
traction lowers the energy, and it follows from the strict submarkovianity that 



the corresponding eigenvalue is strictly smaller than 1. Hence there is a "mass 
gap": For some positive e, the "energy" < ~Lz,z > m dominates e < z,z > m 
for all z. 

2.2 A twisted Gaussian measure 

Then, although L is not symmetric, an elementary computation (given in a more 
general context in the following section) shows that for any \ <= denoting 
M x the diagonal matrix with coefficients given by \-> 

j^p^ J {e~ <z7 ' ^e^< Lz ^mx u dy u = det(-M ro L + M xm )~ l 

. As a consequence, differentiating with respect to Xx, 

7^4m / ^(e" <z "' x>m e5< L ^">™ndx„dy„ = det(-M m L+M xm )- 1 — {-L+M x )~l 

In a similar way by perturbation using a non diagonal matrix, one obtains 
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4m [ z x z^(e- <zz >x > ™e 1 i <Lz > z> ™Ildx u dy u = det^M^+M^y 1 — (-L+M x )^ 



But with the usual notations for Markov processes, setting If — J* A l{ Xs = x } ^—ds 



and I* — l x , we have 



— (-L + M x )~y = E x ( f e~<^>-dl\) 
m y J a 

Defining the path measure E Xi2/ by: E X (J^ G(x s , s < t)dlf) = ~& X , V {G) the 
above relation writes 

— (-L + M x )^=E x , y (e-<^>-) 
m y 

It follows that we have proved the following: 
Proposition 1 For any continuous function F on 

(*) J z x z y F(z u z u ,u e X)e^ <Lz;s>m Udx u dy u 

= J E x ,y{F(l u + z u z u , u e X))e^ <Lz ' z> Tldx u dy u 
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2.3 Positivity 

It should be noted that setting p u = \z^z u and z x = \J p x /2e l6 *the image on R^ 
of the normalized complex measure vx — n^pc\ det(—M m L)e^ <Lz > z>m Ildx u dy u 
by the map z u — > p u is an infinitely divisible probability distribution Q on 
Rf with density * |x| det(-M m L) / e < L sfP* u > e " 9> '»ndft u . Note that the 
positivity is not a priori obvious when L is not m-symmetric. This impor- 
tant fact follows easily by considering the moment generating function ^(s) = 
dct(-L+M ) = det(J + (— L)~ 1 AI S )~ 1 defined for all s with non negative coordi- 
nates, positive and analytic. The expansion in power series around any s (which 
appears for example in [8]) is explicit: 

= detil+i-L+AQ-Hh)- 1 = exp(-log(det(/+(-i+A/ s )- 1 Af, i ))) 
= exp(E t^Tr([(-L + M s )- l M h ] k ). 

As (-L + M s ) _1 is nonnegative, it implies that $ is completely monotone 
as in this last expression, all coefficients of /i-monomials of order n are of the 
sign of (-1)". 

Note that the same argument works for fractional powers of 3>(s) which 
shows the infinite divisibility. Let us incidentally mention it has been known 
for a long time (cf [8] ) that this expansion can be simplified further in terms of 
permanents. 

For x = y, the above proposition then yields the following: 
Corollary 2 For any continuous function F on 

(**) / PxF( Pu , u S X)Q(dp) = J M x , x (F(r + Pu ,ue X))Q{dp) 

Note that this last formula is also obtained in [3] after a direct definition of 
the measure Q. 

Another interpretation of this positivity and of infinite divisibility can be 
given in terms of a Poisson process of loops. It will be developed in a forthcoming 
paper but let us simply mention that Q appears to be equal to the distribution 
of the occupation field associated with the Poisson process of loops canonically 
defined by the Markov chain. 

REMARK: HYcX, it is well known that the trace of the process on Y is 
a Markov process the potential of which is the restriction of V to Y x Y. The 
distribution vx induces vy. Therefore the formulas (*) and (**) on X and Y 
are consistent. 

EXAMPLE: let us consider, as an example, the case where X — {1, 2....N}, 
Qi = 1, ILj = l i<N ,j=i+i, A*i = l;=i, mi = 1 and (-L)rj = 1^.. 



The characteristic polynomial of |(II + II) is (— A + VA 2 — l) N + (— A 

IT 1 



VA 2 — 1)^ hence one gets easily that the mass gap equals 2 sin 2 (27^) 



1 Completely monotone functions in several variables were already used in [T]. 
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Under E X)a all local times vanish except l x which follows an exponential 
distribution. Moreover an easy calculation shows that Q reduces to a product 
of exponential distributions. The formula (**) reduces to the convolution of two 
exponentials. 

3 The infinite case 

We now explain how in certain situations, the above can be extended to a 
Markov process on an infinite space X. Of course, a Markov process for which 
points are not polar can always be viewed elementarily as a consistent system of 
processes on finite subspaces but we aim at a stronger representation allowing 
to consider any functional of the occupation field. There are some obvious 
obstructions to a generalization. The mass gap property does not always hold: 
consider for example the case of a constant drift on an interval, analogous to the 
above example. Some conditions have to be assumed in order that the energy 
controls the skew-symmetric part of the generator. 

3.1 Some calculations in Gaussian space 

Let if be a real Hilbert space with scalar product <, >. At first the reader may 
suppose it finite dimensional and then check that the assumptions we will make 
allow to extend the results to the infinite dimensional case. 

Let <f> be the canonical Gaussian field indexed by H. Given any ONB tk of 
H, Wk = 4>(ek) are independent normal variables. Recall that for all / € H, 

<t>(f) = E fe < f,e k > w k and E(e^) = e"^. 

In the following, </>(/) can be denoted by < <f), f > though of course cf> does 
not belong to H in general. 

Let K be any Hilbert-Schmidt operator on H . Note that K(f> = J2k w kKe k 
is well defined as a i/- valued random variable, and that -E(|| AT</>|| 2 ) = Tr(KK*). 

Let C be a symmetric non negative trace-class linear operator on H. Recall 
that the positive integrable random variable < C<fi, cj) >£ L 1 can be defined by 
J2k < Cek,ek > wf for any ONB diagonalizing C. 

Moreover, £( e -|«^>+#(/)) = det(/ + C)-^e- <<J+C 2 lfJ> (the determi- 
nant can be defined as + K), where the are the eigenvalues of C. 

In fact det(I -I- T) is well defined for any trace class operator T as 1 + 
S^Li Tr(T An ) (Cf [7] Chapter 3). It extends continuously the determinant 
defined with finite ranks operators and it verifies the identity: 

det(7 + Ti + T 2 + TjT 2 ) = det(7 + Ti) det(7 + T 2 ). By Lidskii's theorem, it 
is also given by the product + A^) defined by the eigenvalues of the trace 
class (hence compact) operator T, counted with their algebraic multiplicity. 

Let </>i and <f>2 be two independent copies of the canonical Gaussian process 
indexed by H. Let B be a skew-symmetric Hilbert-Schmidt operator on H . 
Note that < Bcf) 1 ,cf)2 >= < Bek, e; > w\wf = — < Bfo^i > is well defined 
in L 2 and that E(e t<B ^^>) = ^(e-sll^ill 2 ) = det(J + BB*)~? . 
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As B is Hilbert-Schmidt, BB* is trace-class. The renormalized determinant 
det 2 (7 + 7?) = dct((7 + 7?)e~ s ) is well defined (Cf [7]), and as the eigenvalues of 
B are purely imaginary and pairwise conjugated, it is strictly positive. Moreover 
since B* = —B, det 2 (J + B) = det 2 (J- B)= det(I + BB*)i 
Finally, it comes that E{e i<B4 ' 1 ' 4 ' 2> ) = det 2 (7 + B)- 1 
More generally, setting ip = fa + i<f>2, 
E( e -i< c - B ^^ > ) = det 2 (7 + C + B)- 1 exp(-TV(C)) 
(Recall that when T is trace class, det 2 (7 + T) exp(Tr(T)) = det(7 + T)). 
Indeed, E(e^^ <( - C ^ B ^^ > ) = £:( e -5<c*i^i>-5< c '*2,02>+i<-B0i,02>) 
= det(7 + C')-^( e -5«(- f +c)- 1 S0i,B0i>+<C0 1 ,0 1 ») ( by integration in fa) 
= det(7 + C)"3 det(7 + C - B(I + C)- l B)-^ 
= det(7- {I + C)- 1 2B{I + C)- 1 B{I + C)- 1 2)- 1 * + C)- 1 
=det((I+(I+C)-iB(I+C)-i)(I-(I+C)-?B(I+C)-i))-i dct^ + C)- 1 
= det 2 (7 + (I + C)-?B(I + Cyiy^ctil + C)- 1 (since det 2 (7 + (I + 
C)~^B{I + C)~i) = det 2 (7 - (I + C)~iB(I + C)~i) by skew symmetry 
as before) 

= det 2 (7 + C + B)- 1 exp(-Tr(C)). 

Note that I + C + B is always invertible, as C + B is a compact operator 
and —1 is not an eigenvalue. 

Let /i and / 2 be two elements of H. Set D{f) =< /, /i > / 2 . For small 
enough e, £( e -f (<(C-b)^>+sWi)^(/2))) 

= det 2 (7 + C + B + sD)- 1 cxp(-TY(C + sD)) 

= det 2 (7 + C + B)- 1 cxp(-TY(C)) det(7 + e(I + C + B)~ 1 D)~ 1 . 

Hence, differentiating both members at e = 0, 

Emfi M/ 2 )e-5«(c-B)^^») = de^^+C+S)- 1 exp(-Tr(C))Tr((/+C+ J B)- 1 J D) 

Therefore 

£(^(/iM/ 2 )e-^' C -^^') x 

-<(J + C + fl) (h),h>. 

If C is only Hilbert-Schmidt, we can consider only the renormalized "Wick 
square" :< Ccj), cf) >:— J2 k < Cek^k > {wl — 1) for any ONB diagonalizing C 

and £( e -f :<C<M>:+#(/)) = de t 2 (7 + C)4 e - < " +c r" > 
The results given above extend immediately as follows: 

£( e -i:<(c-BW,?>i) = det2 (7 + C + B)- 1 
S(V(/i)^(/ 2 )e-^ :< ( c - B ^^ >: )) 



S ( e -i(:<(C-B)^, 1 /» : ) 



=< (I + C + B)- X {f x ),h > 
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3.2 A class of Markov processes in duality 

Let (V a , a > 0) and V a be two Markovian or submarkovian resolvents in duality 
in a space L, 2 (X,B,m) with generators L and L , such that: 

1) Denoting V = V(L) n T>(L), L(V) is dense in L 2 (m) 

2) < —Lf,f > m > £ < f,f > m for some e > and any / e V (i.e. we 
assume the existence of a spectral gap: it can always be obtained by adding a 
negative constant to L. 

Let H be the completion of T> with respect to the energy norm. It is a 
functional space space imbedded in L 2 (X, B,m). Let A be the associated self 
adjoint generator so that H = D(V '—A). On T>, \{L + L) = A. 

The final assumption is crucial to allow the control of the antisymmetric 
part: 

3) B = (-A) -1 -^^ is a Hilbert-Schmidt operator on H. 

Equivalently, (— A)~ ? L 7, L (— A)~^ is an antisymmetric Hilbert Schmidt op- 
erator on L 2 (to) since for any ONB of H, (— A)^ek is an ONB of h 2 (m). 
Note that I — B is bounded and invertible on H and that Vo = (I— B)~ 1 (—A)~ 1 
maps L 2 (to) into H. Indeed, one can see first that on V, A(I — B) = L so that 
on L(V), Vo = (I-BY^-A)- 1 

EXAMPLES 

This applies to the case of the finite space considered above. 
Let us mention other examples: 

1) Diffusion with drift on the circle: X = S 1 , A = - e, L-L = b(6)^g, 
where b is a bounded function on S 1 . 

Indeed, considering the orthonormal basis e lke in h 2 (d9), -^=e lke is an 

orthonormal basis in H = H 1 , and 

E fc |M^)M^ e i[ 2(de) - £ M V) 2 T^ < - 

2) Levy processes on the circle: The Fourier coefficients at + ibk of L should 
verify £ fc (^) 2 <^ 

3.3 An extension of Dynkin's isomorphism 

Assume X is locally compact and separable, and that functions of H are con- 
tinuous. By the Banach-Steinhaus theorem, given any point x e X, there exists 
an element of H, denoted i] x defined by the identity: f(x) =< r/ x ,f >h ■ Note 
that r)x =y~] e/c(x)e/c for any orthonormal basis of H. 

The resolvent V\ is necessarily Fellerian and induces a strong Markov pro- 
cess. Denote by If the local time at x of this Markov process. Let x and 
y be two points of X. Set: < r] x ,r]y >h= J2 e k(x)ek(y) = K(x,y) so that 
A^f{x) = J K(x,y)f(y)m(dy) or A~ x f = J ■q y f{y)m{dy). 

Set V (x,y) =< I + B)~ lr ri x ,riy > H and note that V (x,y) is a kernel for V . 
Indeed, for any f,g e L 2 (m), < V f,g >m,n)=< {I-B)~ 1 (-A)- 1 f,g > L 2 (m) =< 
(I — B)~ 1 (—A)~ 1 f, A~ x g > H = J f(x)g(y)V (x,y)m(dx)m(dy) 
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As a consequence, Vo(x,y) = E x (l^). 

Applying the construction of the section 3-1, we see the kernel K(x, y) is the 
covariance of a Gaussian process (Z x — ip(ri x ) — J2 e k( x ) w k, x G A), for any 
ONB e k of H. 

More generally, for any non-negative finitely supported measure x = Pj^uj 
on A, letting C be the finite rank operator: C = J2i PjVu, -^Vu, ■,, set V x (x, y) =< 
(I-B + C)- 1 T lx ,r ly > H 

In a similar way as above for V , we have V x (x,y) = E x (e~f l * ^ dz ">dl y t ). 
Then, from section 3-1 

E(Z x Z y )e-^< c - B _^>^) = J E( e -|«(C-B)V>>H)) < (I-B+C)' 1 ^,^ > H 

= £7( e -i(<(c-*)*.*>*))^(z,y) 

But (< Cip,ip > H = J ip(r] u )ip(Vu)x(du) = J Z u Z uX (du) and < Bi/j,ip > H =< 

On the other hand, at least formally in general but exactly in the finite 
dimensional case, 

/ k^k Zu Z u m{du) =< {-A)- 1 ^Z,Z> H =Y.w 1 k w* < (-A)- 1 ^e k ,e l > H 
Hence we can denote: < Bip,ip > H by < ^^Z, Z >L 2 (m)- Therefore 
E{Z x Z y e^ <! ^ z ^Hm)-S 

= E <g> E x ( J C e^^^^Me- fit x(dz)-f z u z uX (du) dl v^ 

Let E* denote the non-normalized law of the bridge of duration f from a; 

to y. Set /x Xi!/ = ^E xy dt, so that /U Xj2/ (l) = E x (l^) — V (x,y) and denote P 

the local time at z for the duration of the bridge. 

Then E{Z x Z v e^ <! ^ z;E> ^^- ! z ^ uX{du) ) 

= Eie i<^z,z> L2(m) -fz u z uX (du) j e -Sl'x{dz) d ^ y) 

Finally, we get that for any continuous bounded function F of N non negative 
real coordinates, and any iV-uple of points Uj in A, 

E(Z x Z y e^ <k ^ z - z> ^^>F(Z Uj Z U] ) = E(g>^ x . y (ei <k ^ Z ' Z> ^^F(l^+Z U] Z U] )) 

and the formula finally extends to 
Proposition 3 For any bounded measurable function of a real field on X : 

(*bis) E(Z x Z y e^ < ^ z ' z> ^(^F(ZZ)) = E® ^(e^^ 2 '^ 2 ^) F(l + ZZ)) 

It induces the formula (*) obtained in the finite case if we consider the trace 
of the process on any finite subset. 

We see also that the restriction of the twisted Gaussian measure 

X <- L — l n ~~y 

e 2 2 ' ^ 2 ^)P(dZ) 

to a(Z u Z u ,u e A) is a probability measure Q under which the distribution of 
the process (Z u Z Ul u € A) is infinitely divisible. Moreover, the important fact 
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is that this probability is absolutely continuous with respect to the restriction of 
P to a(Z u Z u , u € X). It is clear that that formula (**) of the corollary extends 
in the same way to yield the following 

Corollary 4 For any bounded measurable function of a nonnegative field on 
X : 

(* * bis) J p x F( Pu ,u e X)Q{d P ) = J fi x , x (F{l u + Pu ,ue X))Q(dp) 

Hence it follows for example that the continuity of the Gaussian field Z 
implies the continuity of the local time field I under all loop measures p, xx . 

Note finally that these results, as in the symmetric case, can be extended to 
some situations where the local time does not exist (like the two dimensional 
Brownian motion), by considering the centered occupation field and the "Wick 
square" : Z X Z X : (formally given by Z X Z X — K(x,x)) as generalized random 
fields. This makes sense for the Wick square provided if is a Hilbert-Schmidt 
operator. 
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